Electronic transport on carbon nanotube networks: a multiscale computational 

approach 

Luiz F. C. Pereira ^, M. S. Ferreira 

School of Physics, Trinity College Dublin, Dublin 2, Ireland. 



o 

D 
(N 



I 



I 

o 
o 



> 
as 
i> 
m 
in 

o 

m 



X 
J3 



■Abstract 

Carbon nanotube networks are one of the candidate materials to function as malleable, transparent, conducting films, 
with the technologically promising application of being used as flexible electronic displays. Nanotubes disorderly dis- 
tributed in a film offers many possible paths for charge carriers to travel across the entire system, but the theoretical 
description of how this charge transport occurs is rather challenging for involving a combination of intrinsic nanotube 
properties with network morphology aspects. Here we attempt to describe the transport properties of such films in two 
different length scales. Firstly, from a purely macroscopic point of view we carry out a geometrical analysis that shows 
'how the network connectivity depends on the nanotube concentration and on their respective aspect ratio. Once this 
is done, we are able to calculate the resistivity of a heavily disordered networked film. Comparison with experiment 
offers us a way to infer about the junction resistance between neighbouring nanotubes. Furthermore, in order to guide 
,the frantic search for high-conductivity films of nanotube networks, we turn to the microscopic scale where we have 
developed a computationally efficient way for calculating the ballistic transport across these networks. While the bal- 
listic transport is probably not capable of describing the observed transport properties of these films, it is undoubtedly 
■useful in establishing an upper value for their conductivity. This can serve as a guideline in how much room there is for 
improving the conductivity of such networks. 

Keywords: carbon nanotube networks. Green functions, electronic transport 



1. Introduction 

Carbon nanotubes present remarkable physical proper- 
'ties and since their discovery, have shown much promise in 
•terms of applications [1| . Nanotube network films are be- 
ing heavily investigated as the basic building material for a 
new generation of fiexible electronic devices. In particular, 
thin metallic films are currently the most promising mate- 
rial to enable the construction of large scale flexible elec- 
tronic displays Carbon nanotube (CNT) films are 
produced by deposition of solutions containing dispersed 
"CNT bundles. Figure [T] shows a typical nanotube network 
film produced by this technique. In general, it is possi- 
ble to control the density of the films by controlling the 
amount of solution deposited. 

The fabrication of novel electronic displays requires 
■flexible, transparent, highly conductive films and nanotube 
networks seem to possess all the right ingredients. Flexi- 
bility is a basic property of CNTs. Transparency can be 
achieved by controlling the density of nanotubes in the net- 
work. Finally, while individual CNTs are known to have 
remarkable transport properties, the electronic transport 
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in network films takes place through individual nanotubes 
interconnected by junctions. In this case, electrons experi- 
ence low resistance when travelling along a nanotube, but 
are faced with significant resistance at the junctions. Fur- 
thermore, transparency and conductivity are intimately 
connected in thin metallic films. A very dense film might 
show high electrical conductivity but low transparency. 
On the other hand, a sparse network may show high trans- 
parency but poor conductivity. It is essential to search for 
optimum ways of increasing the conductivity but main- 
taining a good level of optical transparency. To achieve 
that, it is necessary to obtain a clear understanding of how 
the electronic transport occurs in such networked films, a 
task in which computer models are extremely useful. 

Modelling complex structures such as carbon nanotube 
network films is a huge challenge. A successful theoretical 
model needs to account for the topology and the disordered 
nature of the arrays formed by individual nanotubes (and 
bundles). It is also necessary to account for the electronic 
coupling between individual components throughout the 
films, since electronic transport takes place through a tun- 
nelling network spanning the whole film. Finally, it is nec- 
essary to account for the electronic structure of individual 
carbon nanotubes. 

The problem presents two clear length scales. One 
which accounts for the geometrical aspects of nanotube 
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Figure 1: Scanning electron micrograph of a typical carbon nanotube 
film produced by deposition of a solution containing dispersed nan- 
otube bundles. Image is courtesy of J. N. Coleman's research group, 
Trinity College Dublin. 



network films, as well as the complex topology of these 
structures. We shall refer to this as the macroscopic scale. 
The other length scale considers that carbon nanotubes 
are made of individual atoms arranged in the usual hexag- 
onal structure. This shall be referred to as the micro- 
scopic scale. Once the relevant length scales have been 
clearly identified, it is possible to develop a multiscale 
model which accounts for both with the appropriate ac- 
curacy. 

The geometrical aspects of individual carbon nanotubes 
can be properly captured by a simple model where these 
are represented by rods of specified length and diameter. 
Choosing a random position and orientation for each rod, 
within a containing volume, a network of interconnected el- 
ements is generated, which mimics the topology of the net- 
works in CNT films. Naturally, percolation theory meth- 
ods can be used to investigate the transport character- 
istics of these networks. On the microscopic scale, the 
electronic structure of a carbon nanotube can be satisfac- 
torily described by simple model Hamiltonians. The tight 
binding approximation, considering only one orbital per 
carbon atom, is enough to provide a description of the 
relevant electronic states. Furthermore, standard Green 
function methods can be used to describe the electronic 
structure and electronic transport characteristics across 
the nanotube networks. 

From a theoretical viewpoint, the description of the 
transport properties of such networks requires a multi- 
disciplinary effort involving classical electromagnetism, quan- 
tum physics and statistics. In this work we attempt to 
use tools of all these multidisciplinary fields of study to 
present a multiscale approach that can describe the elec- 
tronic transport on carbon nanotube films. It begins with 
a description of the model utilised to account for the com- 
plex topological network formed by CNTs. It will be shown 
how the connectivity of the network depends on the den- 
sity of the network as well as the dimensions of individual 
nanotubes. Considering a limit where the network resis- 
tance is exclusively due to nanotube junctions, it will be 
shown how the proposed model compares with standard 



percolation theory predictions. Moving on to the micro- 
scopic description of the problem, it will be shown how 
Green function methods can be used to calculate electronic 
transport on disordered nanotube networks. Finally, by 
combining both length scales in a single formalism, we 
show how the resulting multiscale model can be applied to 
calculate the maximum conductance of network films. 

2. Random Networks of Finite Rods 

The complex structure presented by carbon nanotube 
network films can be modelled by standard geometric per- 
colation models. Single nanotubes and small diameter 
bundles are modelled by rigid rods of specific length and 
diameter. The structure of the film itself can be repre- 
sented by a 3D box which encloses the disordered array of 
interconnected rods. It is important to define a box with 
appropriate dimensions, such that most rods are located 
inside the box. It is also important to reproduce other 
aspects present on the real experiments, such as the place- 
ment of metallic electrodes on the sides of the films for 
electrical resistance measurements. 

Once the aspect ratio of the rods, defined as length- 
to-diameter ratio, and the dimensions of the enclosing box 
are defined, it is also necessary to specify the total number 
of rods inside the box, i.e. the number of nanotubes in the 
film. The total number of rods considered along with the 
dimensions of the box determine the density of the film. 
The density of a specific CNT film can also be expressed 
in terms of its volumetric fraction occupied by nanotubes. 
The volume fraction of a typical laboratory produced nan- 
otube film is related to the volume of CNT solution used 
for deposition Q , and can be directly measured. In a 
computer simulation, the volume fraction depends on the 
total number of rods distributed inside the box, and is 
calculated directly for each configuration generated. 

Each rod is randomly placed and randomly oriented 
inside the box. The arbitrary orientation direction is cho- 
sen in a 3D space with no constraints, and individual rods 
are allowed to overlap with other rods. An arbitrary rod 
might also have part of its volume located outside the box, 
and this portion is not included when calculating the vol- 
ume fraction of the film. Each ensemble of randomly dis- 
tributed/oriented rigid rods gives rise to a complex random 
network of interconnected 1-dimensional elements, capable 
of capturing the most important morphological features 
found on a nanotube network. Figure [5] shows a typical 
configuration generated with N — 300 identical rigid rods 
of length-to-diameter ratio equal to 20. 

It is possible to identify interconnections between rods 
throughout the network by considering the centre-to-centre 
distance between each pair of rods. If the smallest centre- 
to-centre distance between a given pair is smaller than the 
sum of their radii, then this pair is said to be connected. 
By calculating the distance between each possible pair of 
rods in the network, all the connections are mapped and a 
list of neighbours for each rod is constructed. In general. 
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Figure 2: Illustrative representation of rigid rods randomly placed 
inside a rectangular box. There is a total of N = 300 identical rods 
inside the box, each with aspect ratio H/D = 20. 



each rod has an arbitrary number of neighbours, but the 
number of connections (neighbours) per rod follows some 
probability distribution with a well defined mean. The 
average number of connections per node is defined as the 
connectivity of the network. Networks with a low den- 
sity of rods, i.e. a small volume fraction, are expected to 
have low connectivities, while networks with a large vol- 
ume fraction will have high connectivities. In other words, 
the density of junctions in a network film depends directly 
on its volume fraction. 

In order to investigate electronic transport across ran- 
dom nanotube networks, it is necessary to include the 
placement of metallic electrodes in model films. Contacts 
to electrodes are crucial to electronic transport studies, 
and must be addressed with caution. In experimental re- 
alisations, metallic electrodes are placed at opposing ends 
of a film and the resistance is measured between these 
electrodes in a usual two probe method. In this work elec- 
trodes are modelled by two opposing faces of the 3D box. 
A rod crossing one of the specified electrode faces is con- 
sidered to be in contact with that particular electrode. In 
this manner, a map of all the rods connected to the elec- 
trodes is constructed. 

2.1. Simulation Details 

Individual carbon nanotubes are represented by rigid 
rods of length £ and diameter D, with a well defined aspect 
ratio given by £/D. The film is mimicked by a rectangu- 
lar cuboid, whose dimensions are defined in terms of the 
rod length £. In all the results reported here the box has 
dimensions given by 2£ x 2£ x A£. For each rod, a random 
position for one of its ends is generated inside the box, and 
the position of the other end is determined by the orien- 
tation of the rod. The total number of rods considered 
in each film is N, and along with i and D these are the 
parameters of the simulations. 

The orientation of the rods is chosen at random from 
all the possible directions in a 3D space. A particular rod 
can have one of its endings located outside the box, and 
the portion of the rod located outside of the film is not 



considered when calculating the volume fraction. Rods 
are allowed to overlap with other rods, which might lead 
to an overestimation of the volume fraction, especially in 
the case of very dense networks. However, this possible 
overestimation of the volume fraction does not affect the 
morphology of the resulting networks. 

Electrodes are located on the 2£ x 2£ faces of the box, 
and rods crossing one of those faces are connected to that 
electrode. The box is defined with a separation large 
enough to avoid an individual rod being connected to both 
electrodes, in order to avoid short circuits. For a fixed 
number of rods, and a box of fixed dimensions, a very 
large number of configurations can be generated by ran- 
domly distributing the rods. Therefore, a large number 
of independent configurations must be considered, and a 
careful statistical analysis of the results must be carried 
out. All the results reported in the next section are ob- 
tained by averaging a large number of independent reali- 
sations, without discarding any configuration. 

In order to determine which pairs of rods are intercon- 
nected it is necessary to calculate the minimum distance 
between each rod and all the rods around it. Identifying 
which rods are in contact with the electrodes also requires 
testing the criteria for all rods close to the respective faces. 
Both procedures can be very time consuming from the 
computational point of view, and the computation time 
grows very rapidly with the total number of rods consid- 
ered. In all simulations we have considered networks with 
a maximum of 1800 individual rods. 

2.2. Results 

Let us begin by analysing the relation between the to- 
tal number of rods and the volume fraction of the films 
simulated. As previously mentioned, experimental studies 
have found a direct relationship between the volume of so- 
lution deposited and the volume fraction of the resulting 
network 5;]. Figure [3] shows the relationship between the 
volume fraction and the total number of rods randomly 
distributed in the film. The data is shown for different 
aspect ratios, which is controlled by varying the diame- 
ter of the rods. Each data point is an average over sev- 
eral independent configurations and error bars are smaller 
than symbol sizes. Solid lines shown in the graph are lin- 
ear least-square fits to each data set, and confirm a direct 
proportionality between the volume fraction and the total 
number of rods in the network. The existence of such a 
direct relation between the number of rods and the volume 
fraction provides a direct way of comparing and combin- 
ing simulation results with experimental measurements. 
On a real CNT film there are several thousand individual 
nanotubes, whereas computational resources limit simula- 
tions to only a few thousand rods. However, morphologi- 
cal properties such as connectivity and density of junctions 
are directly related to the volume fraction, which connects 
experiment and simulation at least on a qualitative level. 

Now we shall proceed to analyse how the connectivity 
of the network depends on the concentration of rods. For 
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Figure 3: Dependence of the volume fraction on the total number 
of rods. The lines shown in the graph are linear least-square fits to 
each data set, and error bars are smaller than symbol sizes. 
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Figure 4: Average network connectivity as a function of volume frac- 
tion. Inset: Universal scaling of the connectivity with Vf X l/D. 
Error bars are smaller than symbol sizes. 



a film of fixed dimensions and containing rods of identical 
aspect ratio, it is natural to expect that the average con- 
nectivity will increase with the volume fraction. If there 
are more rods inside a fixed finite volume, more rods will 
be in contact. On the main panel of figure 0] it is shown 
how the average connectivity per rod (a), increases with 
the volume fraction V/, for each aspect ratio considered. 
From the data it is also clear that for a fixed volume frac- 
tion, networks composed of high aspect ratio rods have a 
higher connectivity when compared to lower aspect ratio 
rods. This is an interesting feature, which is in accordance 
with experimental observations that films made of high 
aspect ratio bundles present a lower overall resistivity [6| . 

Based on a rod percolation model similar to the one 
considered here, a purely geometric argument has been 
used to propose that the mean number of junctions per 



Figure 5: Connectivity probability distribution for random networks 
of finite rods with aspect ratio 1/ D = fO, and indicated volume frac- 
tion. Vertical solid lines indicate the calculated average connectivity. 
Each data point is an average over 100 independent configurations. 



nanotube (or bundle) is expected to behave as 



{a)=(ix Vf 



(1) 



where /? is a constant independent of the nanotube (or bun- 
dle) diameter and length. Also, Vf is the volume fraction, 
{£) is the average bundle length, and (D) is the average 
bundle diameter. In order to verify if this hypothesis is 
satisfied by the present model, the average connectivity is 
plotted as a function of the scaled variable Vf x £/D, in 
the inset of figure HI Remarkably, the results are in good 
agreement with the predicted relation, for all the aspect 
ratios considered. Besides providing a direct verification 
of equation ^ , from the linear least-square fit on the in- 
set of figure |4] the value of the constant (3 = 2.14 ± 0.03 
is calculated, which is independent of the other factors in 
the equation. 

It is also instructive to analyse the probability distri- 
bution of connections per rod in the network. Figure [S] 
presents the distribution for the probability that an arbi- 
trary rod will have a connections. The data is obtained 
by considering random configurations with rods of aspect 
ratio £/D = 10. The vertical solid lines indicate the cal- 
culated average connectivity (a) for each volume fraction 
specified. From the data it can be clearly seen that the 
average number of connections per rod increases as the 
concentration of the network is increased, while the prob- 
ability distribution becomes broader. The broadening of 
the distribution indicates that as the density of the net- 
works increases, they actually become less homogeneous. 

Another important factor extracted from the computer 
model considered is the number of rods connected to the 
electrodes. In figure [S] (main panel) it is shown that the 
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Figure 6: Number of rods connected to electrodes as a function of 
the volume fraction. Inset: Universal scaling of the electrode con- 
nectivity with Vj X {IjD)^. 



number of rods crossing the electrode faces increases with 
the volume fraction of the film. Similar to the average 
connectivity of the network, the actual increase of {cxe) 
with Vj depends on the aspect ratio of the rods, albeit 
in this case the dependence is more pronounced. In order 
to be useful when compared to experimental data, it is 
desirable to find a universal scaling dependence for (a^) 
in terms of the parameters V/, I and similar to the one 
in equation ([T|). Following an empirical approach, it has 
been found that the relation 



(2) 



provides an excellent collapse of the data on a single straight 
line, as shown on the inset of figure [51 Moreover, the pro- 
portionality constant is found to be fis — 1-35 ± 0.07. 

In practice, equations ([1]) and provide numerical 
values for the connectivity and for the number of connec- 
tions to the electrodes in a model network film, for any set 
of parameters chosen. Furthermore, the expressions ob- 
tained for (a) and {oce) can be used in combination with 
experimental data in order to estimate physical proper- 
ties of real carbon nanotube films Finally, equations 
(dl and ([2]) can also be coupled with electronic transport 
models to provide calculations for the electrical conductiv- 
ity of network films. 

In the next section we present a methodology to cal- 
culate the resistance of a network composed of Ohmic re- 
sistors. This methodology is later coupled with the re- 
sults presented above, to provide calculations for the re- 
sistance of nanotube films from a purely classical point of 
view. Later on, the results presented above will be cou- 
pled with a quantum transport model, providing a more 
realistic microscopic picture of electronic transport across 
carbon nanotube films. 



3. Resistive Netv^rorks 

In carbon nanotube films, the resistance of the nan- 
otubes itself is negligible when compared to the junction 
resistance. In this case, a nanotube network film resembles 
a networks of perfect conductors connected by Ohmic re- 
sistors. Calculating the resistance between arbitrary nodes 
on a lattice of resistors is a problem that has been inves- 
tigated before. Consider for example a resistive lattice 
where each vertex contains a resistance i?o, i-e. between 
every two adjacent nodes there is an Ohmic resistor i?o. 
Different methodologies have been developed in order to 
tackle this problem. 

Venezian Q introduced a method to calculate the re- 
sistance between two adjacent nodes on an infinite square 
grid of equal resistors based on a superposition of po- 
tentials and currents of two one-terminal configurations. 
Atkinson and van Steenwijk [l^ generalised the original 
method to 3D and higher-dimensions. Cserti and collabo- 
rators developed a methodology based on the application 
of lattice Green functions 111, which provides values 
for the equivalent resistance in infinite networks. Build- 
ing on the method developed by Cserti, it is possible to 
propose yet another approach specifically aimed at deal- 
ing with finite networks of resistors. The methodology 
consists on the direct application of Ohm's and Kirchoff's 
circuit laws to finite resistive networks. When a current is 
injected in an arbitrary node and extracted from another 
arbitrary site, the expression for the equivalent resistance 
is given in terms of the potential difference between these 
sites. 

Consider a resistive network where between every two 
adjacent sites there is a resistor i?o- The current that can 
be extracted from a site denoted by A is given by the net 
current from its nearest neighbours, which can be written 
as 



Va+i) + {Va-Va-i) + ' 



(3) 



where Va is the electrical potential on site A, and A±l 
represent two of the nearest neighbours of site A. It is 
possible to rewrite equation ^ as 



CtA 



(4) 



where 5 represents the nearest neighbours of site A, and 
aA its total number of neighbours. In matrix notation, it 
can be written in a compact form as 



I^-^MV. 



(5) 



Only two elements of I are non-zero, the ones correspond- 
ing to the injection and extraction of current, being one 
positive and one negative. 
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The matrix denoted by M has the general form 



M = 



-1 
-1 





-1 

"3 



V 



(6) 



where ua is again the number of nearest neighbours of 
site A, and Ma.b = — 1 if ^4 and B are neighbours or zero 
otherwise. Notice that for a true random lattice, each site 
has an arbitrary number of neighbours. 

At this point, the matrix equation ([SJ is re-arranged to 
read 



V = RqM-H. 



(7) 



By performing the calculation above, we obtain the value 
of the electrostatic potential in each and every node of 
the network. Moreover, from V one calculates the poten- 
tial difference across every resistor in the network, which 
provides a direct way of calculating the equivalent resis- 
tance between any two arbitrary points of the network. 
Furthermore, by constructing the appropriate matrix M, 
it is possible to calculate the resistance between arbitrary 
nodes of any resistive structure, all that is required is to 
describe the matrix M according to the specific network 
connectivity. 

The resistance between two arbitrary sites A and B is 
given by 



Ra.b — 



Va-Vb 



(8) 



Which can finally be expressed in terms of specific ele- 
ments of the inverse of M as 



Roi [M 



A,B 



Ua + [M 



B,B 
B,A )■ 



(9) 



Therefore, the resistance between two arbitrary points in 
any type of lattice can be calculated as long as we know the 
number of connections, and which sites are interconnected. 
It is not hard to generalise this result for the case where the 
resistance between nearest neighbours is not identical, or 
even for a distribution of different resistances. In order to 
do so, the term in Eq. ([3]) would have to be replaced by 
the specific resistance connecting the corresponding sites, 
and could be absorbed into the definition of M. 

3.1. Application 

The method described above can be directly applied 
to model carbon nanotube films where the resistance be- 
tween the two electrodes is calculated through equation 
([9|). Generating a random configuration of rods and cal- 
culating the connectivities of the arrangement, the matrix 
M is written and the resistance between the electrodes is 
calculated. It is important to recall this approach is valid 
when the resistance along the nanotubes is very small, and 
the resistance of the film is dominated by the junctions. 




Figure 7: Schematic representation of the model used to calculate 
the resistance of carbon nanotube networks. Individual CNTs and 
bundles are represented by rods of finite length. Ohmic resistors 
-Ro are place at each junction. Total resistance across the film is 
calculated as the equivalent resistance between electrodes (shaded 
lateral rectangles). 



A random array of nanotubes (represented by finite- 
length rigid rods) is generated inside a containing box of 
specified dimensions, and the volume fraction is calculated. 
The junctions between rods in close proximity are mapped 
and a list of neighbours of each rod is constructed. Sim- 
ilarly, rods that cross the electrode faces are identified, 
and that specific electrode is added to their list of neigh- 
bours. The inter-tube junction resistances experienced by 
conduction electrons, are represented by Ohmic resistors 
(i?o), placed between every two rods in contact. Finally, 
the total resistance between the electrodes is calculated, 
and for the purpose of analysis, it is converted into con- 
ductance by taking its inverse. Figure [7] illustrates in a 
schematic way the set up used in our calculations. 

As previously mentioned, nanotube films are commonly 
produced by deposition of solutions containing dispersed 
CNTs. The disordered network has a density (and volume 
fraction) which is directly proportional to the volume of 
solution used. A small volume of deposited solution will 
generate networks with a low density of connections and 
few percolating paths for the charge to move from one elec- 
trode to the other. In this case, films with high resistance 
will be produced. On the other hand, if the volume of 
solution used is large, it is expected that the films pro- 
duced will have a low overall resistance, since there will be 
several percolating paths for the conduction electrons to 
move across the network. 

It is more convenient to express the charge transport 
results in terms of the inverse resistance, i.e. conductance, 
defined simply as F = 1/R. Percolation theory predicts 
that for small values of V/, the conductance behaves ap- 
proximately as 



(10) 



where Vc is the critical value of the volume fraction. The 
critical exponent t, depends only on the dimensionality of 
the embedding space. In our model the rods are randomly 
oriented in a 3D space, and for this case, the theory pre- 
dicts t — 1.94 [l3|- Beyond the percolation threshold, i.e. 
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Figure 8: Network conductance as a function of the volume fraction 
near the percolation region. Solid lines are a fit of equation l llOll with 
the appropriate 3D exponent t = 1.94. 



Figure 9: Film conductance increases linearly with density of con- 
nections for large enough densities. Log-log plot shows convergence 
to linear regime. 



for networks with a volume fraction eonsiderably larger 
than the critical value 14 1 the conductance is found to in- 
crease linearly with Vf 

The behaviour of T in the percolation region can be 
analysed in detail in order to verify the agreement with 
the standard theory. In figure [SI we fit equation (ITU|) with 
the expected scaling exponent t = 1.94, to the data near 
the percolation region. The agreement with theory is re- 
markable, and from the fit we also extract the value of the 
critical volume fraction for each aspect ratio considered. 
The calculated values of Vc are inversely proportional to 
(■/ which is an intuitive result and is in agreement with 
previous theoretical predictions IJ, |l6|, |l7[ . 

In section [5] it has been shown how the connectivity of 
a random network of rods scales linearly with the product 
of the volume fraction with the aspect ratio of individual 
rods. It is natural to expect the conductivity of the net- 
work to increase with its connectivity, which is indirectly 
demonstrated by the dependence of the conductance with 
the volume fraction (see figurejS]). Moreover, it is also pos- 
sible to analyse the dependence of the film conductance 
with the density of junctions, which is defined as the total 
number of junctions per unit volume, and can be written 
in terms of the network connectivity as 



Nj = 



{a)N 1 
2 V' 



(11) 



where N is the total number of rods inside the box of 
volume V . 

Figure |9] shows a plot of the network conductance as 
a function of {a)N ~ INjV on a log-log scale, and on a 
linear scale on the inset, for three different aspect ratios. 
Once again our calculations are limited by the computa- 
tional resources required for simulations with the rods of 
large aspect ratio. However, it is possible to infer from the 
log-log plot that, at least for high enough concentrations 
of junctions, the conductance seem to converge to the lin- 
ear dependence represented by the dashed straight line. 



Furthermore, in the inset we show the excellent collapse 
of all the data to the same straight line in a linear scale 
T / Rq^ = m X 2NjV. In combination with experimental 
measurements, the model described above has been used to 
estimate the inter-tube junction resistance The result 
obtained is in broa d ag reement with independent exper- 
imental results 



18l . |19|. Nonetheless, it is important to 



point out that in real CNT films the junction resistance 
can vary widely, depending on the chirality of the nan- 
otubes involved as well as the diameter of the bundles, 
and an exact determination of these contact resistances is 
still an open problem. 

4. Microscopic Modelling 

The approach presented in the previous sections anal- 
yses the system in a macroscopic scale. Individual CNTs 
were modelled as perfectly conducting rigid rods of finite 
aspect ratio, and tunnel barriers at nanotube junctions 
were mimicked by Ohmic resistors. The film resistance 
was calculated as the equivalent resistance between op- 
posing electrodes separated by a disordered network of 
identical resistors. As previously mentioned, the method 
provides results in reasonably good agreement with exper- 
iments and previous percolation model studies. 

It is also possible to approach electronic transport prob- 
lems from an atomistic viewpoint, most notably by means 
of ab-initio density functional theory calculations (DFT) 
coupled with non-equilibrium Green function methods |2Ci 
2l|. However, the computational cost of these calculations 
can be considerably high. Moreover, all current DFT im- 
plementations require at least some degree of periodicity 
in order to handle large systems (« 10'^ atoms), which 
unfortunately is clearly not present in disordered CNT 
networks. Nonetheless, ab-initio methods have been used 
to investigate electrical resistance at single wall nanotube 
(SWNT) junctions 2^, and more recent studies have ex- 
tended these results considering the presence of O2 and N2 
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molecules [23| near the junction. The use of first principles 
methods for obtaining the quantum conductance between 
individual nanotubes is able to shed some light on the 
quality of the junction, but it is unable to reproduce the 
general observed features in network films because it is not 
capable of including disorder effects to the scale required. 
The alternative of carrying out a fully atomistic transport 
calculation within a heavily disordered environment is very 
desirable but it is currently too computationally demand- 
ing. It is therefore necessary to compromise if one wishes 
to combine the two features. 

One possible approach consists of using computation- 
ally inexpensive semi-empirical methods to describe the 
electronic structure of carbon nanotubes. Tight binding 
model Hamiltonians are a very efficient and convenient 
way to describe electronic structure of materials in gen- 
eral, and are well suited to deal with carbon nanotubes. 
By describing the electronic structure of individual nan- 
otubes with a simplistic model Hamiltonian, we can afford 
to account for disorder effects by considering hundreds, if 
not thousands, of CNTs, and several possible configura- 
tions for the structure of the random tunnelling networks 
formed in the films. Furthermore, this approach is par- 
ticularly well suited to treat electronic transport across 
networks on a purely ballistic regime. By neglecting all 
sources of decoherence-inducing scattering, it is possible 
to calculate the conductance of ideal network films. The 
conductance calculated in such idealised scenario repre- 
sents an upper bound for the transport performance of real 
films, and can be used as an indication of how much room 
there is for improving the conductivity of these materials. 

4.1. Tight Binding Hamiltonian and Green Functions 

A general semi-empirical tight binding Hamiltonian can 
be written as 



7^1/ = ^ \iy,i)eo{i^,i\ + |i^,i)7o(i^,j|, 



(12) 



where ji', i) represents an electron on atom i of a wire la- 
belled by J/, eo is the on-site energy, 70 is the hopping inte- 
gral parameter, and the sum (i, j) is over nearest-neighbouring 
atoms. This general form can describe several types of 
nanostructures, depending on the choice of atomic struc- 
ture and on the orbital degrees of freedom represented by 
the states Associated with the Hamiltonian T-Li,, there 
are the retarded(-l-) and advanced(— ) single-particle Green 
functions, which are defined as 



(13) 



where 77 is a small positive imaginary part added to the 
electronic energy E in order to avoid singularities at eigen- 
values of the Hamiltonian. 

Electronic transport calculations typically require the 
presence of particle reservoirs, which represent source and 
drain for the conduction electrons. In the case of transport 



simulations of carbon nanotube networks, semi-infinite nan- 
otubes are employed to mimic the presence of particle 
reservoirs. These semi-infinite CNTs are solely described 
by the GF elements relative to the atoms on its final unit 
cell. In order to calculate surface GF elements of semi- 
infinite nanotubes, we apply efficient recursion methods, 
which are numerically stable and computationally inex- 
pensive. 

4-. 1.1. Network Green Functions 

Following the general method applied on electronic trans- 
port calculations, the system is separated into three parts: 
two leads and a central region 2J, 2^. The major differ- 
ence between our particular approach and the usual one is 
that each one of the three parts of the system are them- 
selves composed of several wires. The leads are represented 
by a number of semi-infinite nanotubes, which act as parti- 
cle reservoirs. Meanwhile, the central region of the system 
consists of the disordered network itself, represented by an 
ensemble of finite-sized CNTs. We begin by constructing 
the Green function matrices of a completely disconnected 
system. The inter-tube junctions are introduced in the 
form of electronic hopping terms between CNTs in the 
network, and between the network and the electrodes. Fi- 
nally, the GF of a fully connected network of nanotubes is 
calculated with the Dyson equation. 

Beginning from a disconnected tripartite system, one 
can write the complete Green function as a super matrix 
of the form 



9i 














9i 



(14) 



where are the GF matrix of the left and right elec- 
trodes, and g^ is the GF of the central scattering re- 
gion. At this stage, the network is an ensemble of non- 
interacting nanotubes. 

In this initial description of the disconnected system, 
each of the matrices in 5^ are in block diagonal form, 
with one block for each individual CNT. The electrodes act 
as charge reservoirs which are represented by semi-infinite 
nanotubes. In other words, whenever a CNT is contacted 
to either electrode we will represent it by a semi-infinite 
object. Therefore, there are {as) diagonal blocks given by 
surface GF of semi-infinite carbon nanotubes. In general, 
the electrode GFs can be written as 



S^{E) 







S^{E) 







S^{E) 



V 



,(15) 



where [E) are advanced (or retarded) Green functions 
of the surface elements of semi-infinite wires, which mimic 
the particle reservoirs. 

Similarly, the GF of the central scattering region con- 
sisting of a disconnected network with N wnanotubes is 



8 



represented by a block diagonal matrix with N blocks, 
and it is written as 



9c (E) 











G^(.E) 



V 



(16) 



where we recall that Q^{E) is the Green function of a 
finite-sized pristine monatomic linear chain or carbon nan- 
otube. In general, each diagonal block Q^{E) is a matrix 
whose dimension equals the total number of atoms on the 
individual wires. 

Individual nanotubes are connected by introducing a 
perturbing potential bridging any two CNTs that should 
be in contact. This inter-tube coupling appears in the form 
of an electronic hopping term 7' in the Hamiltonian of the 
complete system, and is given by 



(17) 



where v and v' label the connecting wires, and j and j' 
identify the atoms in each one. Notice that v and must 
be different while j and j' are allowed to have the same 
value. 

These individual hopping terms can be grouped to- 
gether in the form of a general perturbation potential writ- 
ten as 



Vlc 



V = 



V, 



LC 

V, 



Vcc VcB. 



t 

CR 



(18) 







where once again each term represents a matrix itself. 

On the construction of the perturbing potential we as- 
sume that semi-infinite CNTs in the same electrode do not 
interact with each other. Furthermore, we assume that the 
central scattering region is large enough, such that there 
can be no direct hopping terms between semi-infinite wires 
located on different electrodes. The first assumption can 
be easily relaxed, and it should not have any significant 
effects on the results. The second assumption, however, 
should not be relaxed. The physical meaning of the sec- 
ond constraint, i.e., [Vlr] = [Vrl] = [0], is to avoid short 
circuits on the system, which could manifest as spuriously 
high values for the conductance, and significantly affect 
the results. 

The term Vcc is by far the most significant one, and 
also the one that relies at the heart of our methodology. In 
our model, [Vcc] is a matrix constructed from a random 
network, generated from an ensemble of nodes with a given 
average number of connections per node. Equivalently, a 
random graph can also be defined by its total number of 
sites N, and the total number of connections on the graph 
{a)N/2. The process to construct Vcc is as follows: 

1. Select two unconnected wires at random: Vjv' ; 

2. Select one random atom in each wire to carry on the 
interaction: ; 



3. Make {i^,j\VccW'J')=j' ; 

4. Repeat until total number of connections equals (a) N/2. 

A similar procedure is also applied to construct [Vlc] and 
[Vcfl], where the interaction of the electrodes with the 
finite-sized nanotubes is chosen by a randomised process. 
The inter-tube hopping term, represented as 7' is defined 
in analogy with the usual tight binding hopping integral 
parameter 70, and roughly represents the coupling between 
individual atoms in interconnected CNTs. 

Remember that Equations and dT5|) describe the 
GF of the completely disconnected network. The Green 
function of the interconnected disordered systems can eas- 
ily be calculated with Dyson's equation. 



G^=g^+g^V{i-g^Vy'g^, 
and can be written in matrix form as 



Gf G 



G^ 



'L 
^CL 



Gu 



G 



G 



RC 



CJl 

Gr, 



(19) 



(20) 



Equation [20l is the Green function of the whole intercon- 
nected network, including left- and right-hand side elec- 
trodes, central scattering region, and propagators between 
these three regions. Therefore, G^ provides all necessary 
information with respect to the electronic structure of the 
disordered network, and can be used to calculate the elec- 
tronic conductance of the system with the Kubo formula. 

4.2. Kubo Formula Applied to Networks 

The zero-bias conductance of the system is calculated 
with the Kubo formula, which provides a simple expression 
for the conductance by calculating the net electronic cur- 
rent across a reference plane in the system, often referred 
to as the cleavage plane, which is located between any two 
parts of the system. The current is expressed in terms of 
a few Green function matrix elements of the system, and 
the calculation becomes simple enough. 

In analogy to the standard formalism first introduced 
by Caroli et al. 2J], the system is divided in three parts. 



The location of the cleavage plane is arbitrary and can 
be chosen in the most convenient location for the case in 
question. Here it is placed between the central region and 
the right lead, as illustrated in Fig. [TUl Translating Fig. 
[To] into words, regions L and R contain (a^;) semi-infinite 
nanotubes, and region C contains a disordered array of N 
finite CNTs, forming a random network with an average 
number of connections per tube given by (a) . 

Having placed the cleavage plane between the regions 
C and R, the Kubo formula for the zero-bias conductance 
is written as [H, IIS 



nEp) 



4e' 



— j ^e{Tr[GcVcRGRV^j, 
Vc rG RcVc rG Rc]} , 



(21) 
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Figure 10: Schematic representation of the system, divided in left and 
right leads (L,R), and a central scattering region (C). The cleavage 
plane is placed between C and R. 



Therefore, because of the particular structure in Dyson's 
equation, it is possible to reduce the Green function ma- 
trices, such that only the elements affected by the pertur- 
bation V are stored. Consider, for example, a CNT with 
several thousand atoms. If only 20 of these atoms interact 
with neighbouring nanotubes, then only the elements cor- 
responding to these atoms and the propagators between 
them need to be stored, since only these will be affected 
by the perturbation. In this case we can perform calcu- 
lations up to a thousand wires, each with a few thousand 
atoms. 



where all the GFs are calculated at the Fermi energy Ep, 
and G{E) are defined as 



(22) 



4-2.1. Reducing Green Function Matrices 

Even though it is possible to write an explicit expres- 
sion for the Green function of the interconnected network 
as in Equation (1201) . from which the matrix elements re- 
quired by the Kubo formula in Equation (1211) can be ob- 
tained, it can still be an extremely challenging task. The 
difficulty arrises from the size of the matrices involved. 
Consider for example, a network with N — 100 wires, 
where each wire has L — 1000 atoms. In this case, 
is a matrix of {N x L)^ — 10^" elements, and [Vcc] is a 
matrix of same dimension. Storing matrices of this size on 
the memory of computer clusters becomes a demanding 
computational task in itself, and inverting matrices of this 
dimension takes a considerable time, even if one applies 
advanced linear algebra routines. 

The solution to this limitation is facilitated by two gen- 
eral properties of the system. The matrix g^ is block di- 
agonal and y is a sparse matrix. Therefore, the prod- 
uct g^V on Equation (fT9)) is also sparse, or at least has 
many vanishing terms that will not affect G^, i.e., several 
of the matrix elements in do not depend on the per- 
turbation, whenever the corresponding element vanishes. 
In practical terms, the corresponding row/column in the 
matrix G^ will not involve any ^-dependent terms and 
is simply given by Q^. A connection between atom j of 
tube v and atom j' of tube v' will only involve the corre- 
sponding matrix elements of the (^-matrix. Therefore, if a 
given CNT is connected to three neighbours at intra-tube 
sites i, j and fc, the matrix elements Qf^^^, Ofj, ^^j., G^j , 

Stk ' (s-nd the respective adjoint elements) are the only 
matrix elements appearing in the expression for the net- 
work GF of Eq. p^ . This can be an enormous simplifica- 
tion when compared with the full diagonalisation method. 
Rather than inverting enormous matrices of dimension of 
N X L, this technique allows us to express the network 
GF in terms of inverse operations of much smaller matri- 
ces, and whose sizes are defined by the total number of 
connections in the system, i.e. matrices of order N x (a). 



5. Multiscale Model of CNT Networks 

The results from the geometrical modelling presented 
in section[5]can now be combined with our quantum trans- 
port model, providing a multiscale approach to the study 
of electronic transport in carbon nanotube network films. 
In particular, by considering an idealised ballistic trans- 
port regime it is possible to evaluate the upper limit for 
the conductance of CNT network films. It has been ar- 
gued that the construction of transparent, flexible dis- 
plays requires thin films with an electrical conductivity 
o-dc > 6 X 10^ S/m [11]. However, the production of CNT 
films with such high conductivities has been very difficult 
to achieve, and it is possible that nanotube networks are 
not the most suitable material to be used in the next gen- 
eration of flexible displays 2^, 33| . 

In section [5] results for the connectivity of disordered 
networks of rods were presented. Universal relations for 
the dependence of the network connectivity with experi- 
mentally accessible quantities such as the volume fraction 
of the fllms and the aspect ratio of individual wires, were 
obtained. In particular, it was found that for a random 
network of rods with individual length £ and diameter D, 
inside a rectangular cuboid fllm of dimensions given by 
2£ X 2£ X At, the average connectivity per rod is given by 



(a) = (2.14 ±0.03) XT//-. 



(23) 



Furthermore, by considering that metallic electrodes are 
located on the opposing square faces of the fllm, and that 
any rod crossing these surfaces is taken to be in contact to 
that specific electrode, we found that the number of rods 
crossing each electrode can be written as 



{aE) = (1.35 ±0.07) X V> ( 



(24) 



Therefore, once the aspect ratio of the wires and the vol- 
ume fraction of the film are chosen, (a) and {ole) are auto- 
matically defined. These parameters are then fed into the 
general quantum transport model presented in the pre- 
vious section, which yields the conductance for a network 
with the specified geometrical properties. Finally, since we 
have considered disordered networks to be defined inside 
a film of specified geometry, it is then possible to convert 
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the network conductance into film conductivity. A result 
which can be directly compared with experimental mea- 
surements. 

Let us begin by investigating the conductance of net- 
works made of (5, 5) carbon nanotubes, which have a metal- 
lic band structure. The diameter of a (5, 5) CNT is w 0.68 
nm, while the length depends on the number of unit cells 
considered. In a (5, 5) SWNT, each unit cell has 4 x 5 = 20 
atoms. We consider nanotubes of aspect ratio equal to 
11,16.5,22,27.4, which correspond to lengths of 7.5 nm, 
11.2 nm, 14.9 nm, 18.6 nm, respectively. For each of these 
aspect ratios, it is possible to vary the total number of 
carbon nanotubes in the films, such that the network vol- 
ume fraction ranges from 0.065 up to 0.325. Moreover, 
through equations and the combinations oi £/D 
with Vf provide networks with (a) in the range 3.6 — 7.5, 
and (as) = 27 - 77. 

In figure [TT] it is shown how the network conductance 
depends on the volume fraction for the respective aspect 
ratios considered. Each data point is an average over at 
least 50 independent configurations, and the error bars 
represent the statistical error in the ensemble. Similarly 
to the resistive network results presented in section 13.11 
there is a clear improvement in the conductance with a 
higher concentration of nanotubes. Moreover, it is also 
clear that higher network conductance values are achieved 
when the film is composed of nanotubes with larger aspect 
ratios. The dependence of the film conductance on the 
aspect ratio of individual nanotubes is in agreement with 
previous experimental studies 31, S^]. One possible way 
to account for this trend is to consider that, on average, 
a longer nanotube will have a larger distance between its 
junctions with other tubes. Therefore, on nanotubes with 
higher £/D ratio the electrons will, on average, travel a 
longer distance along the tubes before experiencing heavy 
scattering caused by the inter-tube junctions. 

It is also possible to analyse how the conductance de- 
pends on the total number of junctions in a network film 
of volume V, which is defined as NjV = {a)N/2. Figure 
1121 shows the averaged network conductance as a function 
of the total number of junctions in simulated films. The 
dashed lines are least-square fits to each data set, and show 
a linear dependence of the conductance with the number 
of junctions. We again see an increase in the overall film 
conductance as the aspect ratio of individual nanotubes is 
increased. Moreover, the model reproduces the approxi- 
mate linear trend observed in experimental reports [8|. 

It is possible to convert the calculated network conduc- 
tance into film conductivity by considering the geometry 
of the films. Employing the usual relation between con- 
ductance and conductivity, it is possible to write 
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Figure 11; Conductance versus volume fraction for networks of (5, 5) 
armchair SWNTs. The legends indicate the aspect ratio of each 
data set. Higher aspect ratio CNTs provide better conductance, in 
accordance with experiments l3l| . 
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Figure 12: Averaged film conductance as a function of total number 
of nanotube junctions. Dashed straight lines are least-square linear 
fits to simulation data, in accordance with experimental results of 
ref. g. 



(25) 



At this point, it is important to remember that F rep- 
resents the conductance of an idealised network of finite- 
sized carbon nanotubes, in a fully ballistic electronic trans- 
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port regime. In the calculation of the conductance, all 
decoherence-inducing factors (certainly present in experi- 
mental realisations) were neglected. The calculated con- 
ductance is an upper limit for the conductance of a real 
film made of similar CNTs. The conductivity calculated cr, 
also represents an upper bound for the maximum electrical 
conductivity of real nanotube network films. 

Based on an empirical analysis of the conductivity data, 
it is found that it scales with a combined variable of the 
parameters Vf, i and D. This dependence allows for a 
direct comparison of experimental data against this supe- 
rior limit estimate. A comparison between experimental 
and simulation results, provides a direct indication of how 
much it is possible to further improve the transport perfor- 
mance of the laboratory produced films. Defining a scaled 
variable x = Vf x £/D'^, figure [T51 shows the dependence 
of the film conductivity (in units of Siemens per meter) 
with X. The data points represented by circles correspond 
to nanotubes of varying aspect ratio, and includes all the 
data in figure 1121 The data corresponding to networks of 
finite-sized (5, 5) SWNTs shows considerable dispersion, 
even after averaging over several configurations. Nonethe- 
less, it is possible to identify a general trend on the data 
which stresses the linear dependence of the conductivity 
with X, represented by the least-squares linear fit indicated 
by the dashed (black) line. The dispersion on the data 
for CNTs is caused by the finite-length of the nanotubes 
considered. Since the Green function of independent nan- 
otubes is calculated numerically, it quickly becomes too 
computationally demanding to increase the length of the 
nanotubes much further. 

In order to overcome this computational limitation, one 
can investigate the conductance of networks made of ide- 
alised nanowires, represented by linear monatomic chains. 
In order to compare the results with carbon nanotubes, an 
orbital degeneracy is introduced in the chain Hamiltonian 
such that there are two independent conducting channels. 
These monatomic chains present a maximum conductance 
of Ae'^/h at the Fermi energy, just like carbon nanotubes. 
It is possible to write an analytic expression for the Green 
function of linear chains 33| , and very long chains can be 
simulated with little computational effort. In order to min- 
imise the finite-size effects on the local density of states, 
chains with L — 16000 atoms where considered. 

Also shown in figure[T3]is data obtained for the conduc- 
tivity of networks composed of ideal nanowires (squares), 
with varying volume fraction and aspect ratio. In this 
case we can see a much smaller dispersion on the calcu- 
lated conductivity. Furthermore, by fitting a straight line 
to the nanowire conductivity we find 

an = M X 4.25 X IQ-^ xVf^ (in S/m), (26) 

where the factor M accounts for the number of conduction 
channels in the wires, for the case considered here M = 2. 
Because of all the idealisations applied to our calculations 
of the conductivity, the expression above is deemed to be 
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Figure 13: Maximum conductivity of carbon nanotube films as a 
function of scaled variable x = Vf X £/D^. 



the upper limit for the conductivity of random network 
films. 

It is instructive to test the expression in equation (j26p 
for typical SWNT values, namely £ — 1/im and d ~ 1.2 nm. 
For Vf ~ 30%, the predicted upper bound would be au = 
1.8 X 10''S/m if CNTs of these dimensions could be fully 
dispersed to form the network. However, carbon nan- 
otubes are known to bundle together, which means that 
in reality wire diameters are considerably larger. On the 
other hand, larger-diameter bundles have more nanotubes 
on the surface leading to more current-carrying channels 
per wire. Taking all this into consideration, we can com- 
pare our expression with the highest-conductivity case re- 
ported so far (Vf = 30%, ^ = 5 /im, d = 20 nm) [J, [H Q . 



Our prediction of cr„ = 9 x 10^ S/m is only one order of 
magnitude superior to the measured value of ct = 6 x 
10^ S/m. Bearing in mind that the upper bound here ob- 
tained assumes a number of ideal conditions that are ex- 
perimentally unavoidable, this might be a clear indication 
that we are approaching a saturation point in the conduc- 
tivity of nanotube network films. 

Finally, although our focus has been on disordered net- 
works comprised of carbon nanotubes, we can extend our 
results to deal with other wires. This could represent the 
case of networks made of other conducting materials, such 
as noble-metal nanowires, for instance. In this case the 
number of conducting channels M depends linearly on the 
wire diameter and the overall conductivity of the network 
is likely to scale inversely with £), with a proportionality 
constant that depends on the specifics of the material in 
question. If the conductivity of CNT network films is ap- 
proaching its saturation point, it is likely that wires other 
than nanotubes may occupy the post of ideal components 
for disordered network films. Indeed, it has been shown 
recently [19] that silver nanowire network thin films might 
be a strong contender in the race for the next generation 



12 



of transparent electrodes. 

6. Concluding Remarks 

In this work we have shown how a relatively simple 

methodology can be applied to model the complex features 
present in the morphological structure of carbon nanotube 
network films. A multiscale approach has been applied to 
deal with the despairingly different length scales present. 
From a macroscopic point of view, individual nanotubes 
are modelled as rigid rods of finite length and diameter. A 
network of interconnected rods represents the film and the 
volume fraction of these networks allows for a connection 
between theory and experiment. On the other hand, from 
a microscopic approach carbon nanotubes are described by 
a semi-empirical Hamiltonian, and Green function meth- 
ods can be applied to calculate the electronic conductance. 
Merging both scales in a uniqiie multiscale approach is pos- 
sible when the connectivities derived from the macroscopic 
model are used as input into the microscopic conductance 
calculation. 

The volumetric fraction of the film occupied by rods 
was found to be linearly dependent on the number of rigid 
rods in the network. The average number of connections 
per rod in a random network was found to scale universally 
with a combined variable given by the product of the vol- 
ume fraction with the aspect ratio of the individual rods 
considered. Similarly, the number of rods connected to 
electrodes placed in opposing faces of the film was found 
to scale with the volume fraction times the square of the 
aspect ratio. We have presented a method, based on previ- 
ously developed methodologies, to calculate the equivalent 
resistance between arbitrary points on a network of resis- 
tors. The rigid rods percolation model was combined with 
resistance calculations in order to obtain the conductance 
of model carbon nanotube network films. We have shown 
how our results are in perfect agreement with standard 
percolation theory, and how the conductance; of nanotube 
films is afii'ected by the morphology of the networks. Direct 
comparison of simulation results with experimental mea- 
surements are found to be in agreement, and provide es- 
timates for the junction resistance in real world nanotube 
films. These estimates are compatible with independent 
measurements reported in the literature. 

We have shown that an electronic structure descrip- 
tion based on semi-empirical tight binding Hamiltonians 
is capable of accounting for the heavily disordered envi- 
ronments found on nanotube networks. In conjunction 
with results from the geometrical percolation model, we 
have presented simulation results for the conductance of 
single-walled carbon nanotube networks. Finally, by con- 
sidering a fully ballistic transport regime on films of ideal 
nanowires, and comparing the ensuing results with typi- 
cal experimental evaluations we have shown that carbon 
nanotube network films might be approaching their perfor- 
mance limit, and that perhaps metallic nanowires repre- 
sent a more promising alternative to future developments. 
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